Approximate Solutions for the Coupled 
Line Equations 
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The coupled line equations for only two modes, representing the TEoi 
signal mode and a single spurious mode in circular waveguide, are solved 
in series form by the method of successive approximations. Bounds are 
Sound on the magnitudes of the terms in the series solution. These bounds 
decrease rapidly only for "short" waveguides; for long guides many terms 
of the series must be included in the solution. 

The coupled line equations are transformed to a new form, in which one 
of the unknowns A is given by A = — In Go , where Go is the (complex) 
TEoi transfer function of the original coupled line equations. Thus Re A = 
— In | Go | , the TEoi loss in nepers, Im A = — Z G , the TEoi phase in 
radians. These transformed equations are again solved by successive ap- 
proximations; the first term is the commonly used solution that has been 
obtained by physical arguments. Bounds are determined for the magnitudes 
of the terms in these series solutions; for a suitable restriction on the coupling 
coefficient that includes many cases of practical interest, these bounds 
decrease rapidly for long guides. 

In present calculations of the TEoi loss statistics in random guides, only 
the first term of the series expansion for A is considered. Unfortunately 
this approximation has not so far been justified. 

I. INTRODUCTION 

Consider the coupled line equations: 

U{z) = -r / (2) +M*)ii(*), 

h'{z) = +je(z)h(z) - r,/,( 2 ). 



(1) 



These equations are of interest in many applications. Our particular 
interest in them in a companion paper 1 is that they describe the effects 
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of coupling between the TE i signal mode, represented by the complex 
wave amplitude I , and a single spurious mode, represented by the 
complex wave amplitude h , caused by geometric imperfections in 
circular waveguide. We have, of course, assumed that only a single 
spurious mode has significant magnitude, so that all other spurious 
modes may be neglected. For example, we may consider copper wave- 
guide with a rather general straightness deviation; the most important 
spurious mode under many conditions will be the forward TE12 (both 
polarizations must, of course, be considered unless the straightness 
deviation is confined to a single plane). However, these equations apply 
to a variety of other problems which may be described by only two modes 
with varying degrees of accuracy. 

In copper waveguide if the wall losses may be neglected the propaga- 
tion constants T and I\ are pure imaginary and the coupling coefficient 
c(z) is pure real. In helix guide, where loss is added to the spurious mode, 
the propagation constant Ti has a significant negative real part; further, 
as shown by H. G. linger, 2 the coupling coefficient c(z) also becomes 
complex. 

The case where the geometric imperfection (e.g., straightness devia- 
tion) and hence the coupling coefficient is a stationary random process, 
perhaps Gaussian, is of great interest; here it is desired to compute the 
statistics of the TE i transmission I in terms of the statistics of the 
coupling coefficient c(z). Since exact solutions to (1) are easily found 
in only a few special cases, this has been done by using an approximate 
solution to these equations that is essentially a second-order perturba- 
tion solution, and by studying the statistics of this approximate solu- 
tion. 1 The present paper will discuss this approximate solution, will give 
some bounds on the convergence of the approximation, and will indicate 
a basic gap in our knowledge concerning this problem. 

Equation (1) represents a drastic idealization of the real TE01 trans- 
mission problem, in that it contains only one spurious mode and neglects 
all other spurious modes. The approximate solution includes all second- 
order terms; a physical interpretation of this solution states that con- 
version of TE01 to each spurious mode and subsequent reconversion to 
TE01 is considered at all pairs of elementary mode converters, but that 
higher-order terms involving more than one pair of elementary mode 
converters are neglected. 1 The exact solution of (1) includes all higher- 
order terms involving the single spurious mode, but neglects many more 
higher-order terms involving the many other spurious modes that have 
been neglected in (1 ). In view of this it may appear questionable to try 
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to deal with ( 1 ) in more exact terms for the TE i mode conversion 
problem. However, a start has to be made somewhere, and it seems 
unlikely that the general case involving an infinite number of modes will 
be understood before the two-mode case of (1) is understood. Even 
this simple idealized case does not yet have a really satisfactory treat- 
ment. Also, ( 1 ) does apply more or less rigorously to many other situa- 
tions than the TEm mode conversion problem. 1 

In dealing with these equations it is convenient to introduce the 
following change of variables: 



h(z) = e~ r ° : -G (z). 
Then ( 1 ) becomes: 



h(z) = e^r-Gtiz). 



(2) 



(3a) 



G 9 '(z) =jc(2)e 4r '(?i(*), 

Gx(z) = jc(z)e- AT *G (z). 

AV = Aa + ./A0 = To - T, ; 

Aa = a — «i < 0, (3b) 

A/3 - /3„ - /3, . 

Note that we assume Aa < 0, because in circular waveguide the TE i 
signal mode will have lower heat, loss than any of the spurious modes. 
We will always take as initial conditions at z = a TE i wave of unit 
magnitude and zero phase, and a spurious mode of zero magnitude: 

Go(0) = 1, <n(0) = 0. (4) 

II. SOLUTION OF THE COUPLED LINK EQUATIONS BY SUCCESSIVE APPROXI- 
MATIONS (PICARI)'S METHOD) 

We summarize the solution of (3) by successive approximations. Let 
Gu 0l) (z) and (V|<„)U) be the /?"' approximation to the solution of (3). 
Let the initial approximation be given simply by the initial conditions 
of (4): 

G m) (z) = 1, G m (z) = 0. (5) 

Then following Picard's method' 1 ' 4 we obtain the successive approxima- 
tions as follows: 
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G m) (z) = 1, G im (z) = 



G m) (z) = l+j f c(s)e Ars G 1 (o ) (s) 

JO 



da 



= 1 



OmM =3 [* c(s)e-* r 'G 0(0) (s) 

J o 

= i / c(s), 

Jo 



<lx 



e- aia ds 



Gom(z) = 1 +J / c(«)«X(«) & 
Jo 

= 1 - f c(8)e AT 'ds [° c(t)e- Ari <lt 
Jo Jo 

Gi i2 )(z) = j I c(s)e _Ar 'Gou)(s) rfs 

Jo 
= j I c(s)e~~ A " ds 



(6) 



Go ( „)(«) = 1 +j [ c(s)e Ar 'Gi (n -i)(s) ds 

Jo 

GiwOO =i rc(s)e- AT 'Go i n-i)(8)d8. 

Jo 

The w ,h approximation is obtained by substituting the (n — l) th ap- 
proximation in the right-hand side of (3) and integrating. Writing 



wo have 



G 0{n) (z) - Go in -i){z) = go(n)(z), 
G Hn) (z) - Gi in -i)(z) = gwo(z), 

n 

Gom(z) = 1 + ]&0DCft)(*)i 
fc=l 

Gi(„)(z) = F »!(*)(«), 



(7) 

(8a) 
(8b) 
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where the g's are given as follows: 

ffoooU) =j f c(«)e AI W- U («) <fo, « ^ 1. (9a) 

0«.>(*) =j [ c(s)e- Ar Vo (n _i)(s) ds, n 2: 1. (9b) 

•'o 

0o<o)(*) = 1, <7i(o,(z) = 0. (9c) 

It is readily seen that 

0ooo(z) = 0, n odd 

0«n)(z) = 0, neven 



(10) 



so that only even terms appear in the summation of (8a) and only odd 
terms in the summation of (8b). 

Id the standard proof of Picard's method the series of (8) are shown 
to converge to the unique solution of the coupled line equations, (3), 
and bounds are given on the magnitudes of the terms in (8). Thus we 
may write 

Oo(z) = 7L9om(z), (11a) 

»=0 
oo 

(h{z) = E&Hn)(2), (Hb) 

h = 

where the g'a are given in (9) and (10). However, better bounds than 
those given by Picard's general method may be found for the present 
special case. We show that 






Jo 


c(s) 


ds\ 






n\ 













, n odd. 









, n even. 



0K»)(z) 



f | c(s) I ds] 
^ U^ 1 e - Sa: , n odd. 



(12b) 



Suppose that ( 12a) is true for some even value of n. Then from (9b) 
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If \c(s)\dsJ 

Aat \_ J J 



| QHu+M I ^ / I c(0 

n\ Jo 



dt 



o n\ 

c(s) \ds\ d\ J \c(s) \ds\ (13) 






r/ o * i c( s ) i dsi 



— Aaz 



(n+1)! 
where we recall from (3b) that Aa < 0. Substituting (13) into (9a), 



| <7o<„+2)(*) I ^ [ I C(0 

Jo 

^ 1 

" (» + 



Aa( -in! 



/ I c(s) I ds 
Jo J dt 



n + 1)! 

In+l 



hi^[0 c(s)ld8 T d [0 cis)ld8 ] (14) 



£ I c(s) I dsj 



(» + 2)I 

Noting (9c), the results of (12) hold for all n by induction. 

We may ask whether the bounds of (12) are the best that can be 
obtained in general, or if by being sufficiently clever we can do better. 
It is easy to find examples whose terms are actually as large as those 
given in (12), so that no improvement in these results is to be expected 
unless suitable restrictions are placed on the problem. Thus, consider 
the following special case : 

Ar = 0. (15) 

The coupling coefficient is non-negative but otherwise arbitrary. The 
general solution to (3a), subject of course to the initial conditions of 

(4), is 1 



Go(z) = cos 



/ c(s) ds . 



Expanding the cosine in power series, 

J c(s) del \ f c(s) ds 



GM =1-4= 



+ 



(10) 



(17) 



2! " 4! 

The successive terms of (17) are simply the gaM given in (9) and (10). 
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It is readily seen that the magnitudes of these terms are equal to the 
bounds given in (12a) if we require ^ c(z) so that | c(z) | = c(z). 

As another similar example, let AT t± and c(z) be a single 5-function 
located at Zq , 

c(z) = C-Uz - 2o). (18) 

In our present case, i.e., straightness deviation, a discrete coupling of 
the type given in (18) corresponds to a discrete tilt located at z = z a . 
The solution to (3a), subject again to the initial conditions of (4), is 

G (z) = cosC, z > 2„. (19) 

Expanding the cosine, 

<M.) -!-£ + £--. (20) 

Again the terms of (20) are the g (_ n ) of (9) and (10), and their magni- 
tudes are equal to the bounds given in (12a). Of course, this above solu- 
tion, which mathematically is valid for an arbitrarily large tilt in the 
present idealized two-mode case, must fail for large tilts in the physical 
case, the error being caused by neglecting the higher-order spurious 
modes excited by the tilt. While this serves as a further warning against 
uncritical application of the results of the two-mode theory to the 
physical problem, it is still of interest to inquire into the mathematical 
properties of the solutions to (3). 

It is often desirable to express the TE i loss in db rather than as the 
magnitude of the TE i normalized gain, \G a \ . Define the complex TE i 
loss A as 

A = -In Go = A -jO. (21) 

Then 

A = -In | Go | , 

(22) 

e = z Go . 

A is the TEni loss in nepers; the TE i loss in db is simply 8.686 A. If 
we have a number of sections of waveguide separated by ideal mode 
filters, the over-all TEoi gain, Got , and loss A T , are given by 

Got — «ro,i "0,2 ' ' ' j 

A T = Ax + A, + ■ • • , 

(23) 
A T = A, + A % + ■•• , 

e T = e, + e 2 + • • • . 
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The statistics of A and 9 for the composite guide may thus be expressed 
simply in terms of the statistics of A and 9 for the individual guide 
sections. 

Suppose that the series solution for G (z) given in (11a) converges 
very rapidly, so that only the first two terms need be retained. Then 
we have from (8-11) 

G (z) pa G m (z) (24) 

so that approximately 

G (z) = 1 - f c(s)c Ar *ds f c(t)e- Art dt. (25) 

•"o *'o 

Then, assuming that the second term is small compared to 1, we have 
approximately : 

| G (z) | = Re Go = 1 - Re ]f (26a) 

A = -In Go = JJ (26b) 

A = -In | Go | = Re JJ (26c) 

9 = Z Go = -Im// (26d) 

where // is shorthand for 

ff = j 1 c (s)e Ar * ds f c(t) e~ Ari dt (27a) 

= f" e Ar " du [ " c(s)c(s + u) ds (27b) 

Jo •'0 

= -('[' c(s)c(0c Ar """ ds dt. (27c) 

2 Jo Jo 

If the coupling coefficient c(z) is pure real but the differential propaga- 
tion constant is complex (possibly not a physical case), (26c) becomes, 
using (27c), 

A = -In I Go I = f e Aau cos A/3m du f c(s)c(s + u) ds. (28) 
Jo J o 

If c(z) is complex, it turns out that for uniform waveguides its phase 
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angle remains constant and only its magnitude varies, so that we may 

write" 

c(z) = (Cr+j <?,)«(*), (29) 

where c(z) is real. Then (28) becomes 

A = -In [ Go | = (C r 2 - &*) f ('*"" cosAfiudu 

• / c(s)c(s + u) ds - 2C r Ci / e Sa " sin A0u du (30) 
Jo J u 

rL-U 

• I c(s)c(s + «) ^s. 

Jo 

These approximate expressions of (20-30) may be regarded as the 
first terms of series expansions for the various quantities. The above 

approximations will be valid when // « 1 and when the higher- 
order terms may be neglected. From the above analysis it would appear 



that when 



» 1 , all of the above approximations would fail, since, 



in particular, (26a) obviously fails. In spite of this fact, (26b-d), (28) 
and (30) may remain valid for a wide class of long guides of practical 
interest; roughly speaking, the required conditions are that the differ- 
ential loss [ Aa | be large enough and that the coupling coefficient c(z) 
be sufficiently small and uniformly distributed in an appropriate sense. 
This result has been suggested by simple physical arguments; a formal 
mathematical derivation starting with the appropriate restriction on 
c(z) and Aa is given in the following section. These results are of im- 
portance because in a random guide the expected value of the / / term 

increases linearly with distance z\ while the approximations of (25) 
and (26a) fail, the results of (2Gb-d), (28) and (30) may remain valid, 
and so provide us with a theory for long guides. 

It is apparent that further restrictions are required to obtain these 
additional results, by considering the example of (18-20). Thus, let 
the magnitude of the 6-function coupling coefficient be tt/2, so that we 
have in (18.) 

C = T - t c(z) = l-hiz - z ). (31) 

Then from (19) 

G (z) = 0, z > z (32) 
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so that 

A = -In G = oo. (33) 

However, the approximate result of (26b) yields 

A = -InGo^ = jp (34) 

The approximation of (34) is obviously invalid; since these relations 
are independent of Aa, this approximation remains invalid no matter 
how high | Aa | becomes. Cases such as this are ruled out by the addi- 
tional restrictions that require the coupling coefficient to be more or 
less uniformly distributed with z in a certain sense, described in the 
following section. 

III. TRANSFORMATION OF THE COUPLED LINE EQUATIONS TO LOGARITHMIC 
FORM, AND SOLUTION BY SUCCESSIVE APPROXIMATIONS 

We repeat for convenience the coupled line equations, given in (3), 
together with the desired initial conditions, (4). 

G ' (z) ^jc(z)e Arz G 1 (z), 

G.'iz) =jc(z)e-* Tl Go(z). 



(35) 



G (0) = 1, d(0) = 0. (36) 
Next, the following transformation of variables is made : 

G (z) = e~ M2) . (37a) 

Gy(z) = e~ Uz) -H{z). (37b) 

The transformation of (37a) is dictated by the desire to obtain a series 
solution for A, denned in (21). That of (37b) was obtained partly by 
trial and error and partly by intuitive means. Substituting (37) into 
(35), we obtain: 

A'( 2 ) = -jc(z)e ATz H(z) (38a) 

H'(z) = j c(z) e ATz + A'(z) H(z). (38b) 

By substituting (38a) into the second term on the right-hand side of 
(38b), we have: 

A'(z) = -j c (z) e AVz H{z) (39a) 

H'(z) = j c(z) e- AT ' -j c(z) e AVz H\z). (39b) 

The initial conditions of (36) transform via (37) to 

A(0) = 0, #(0) = 0. (40) 



APPROXIMATE SOLUTIONS 1021 

The method of successive approximations may now be applied to (39) 
(or 38), subject to the initial conditions of (40). We note that A(z) is 
absent from (39b), so that this equation contains only a single dependent 
variable, H(z). Thus the successive approximations to H(z) may be 
found without reference to (39a) or to \(z) ; the corresponding approxi- 
mations for A(z) are then found by a simple integration of (39a). We 
note further that (39b) for H{z) is a Riccati equation. 

Thus, let A„(z) and H„{z) be the n th approximation to the solution 
of (39), subject to the initial conditions of (40). Then: 

ffo(«) = 

H } (z) = j f c(s)e" Ar ' ds 

H 2 (z) = j J c(s)e~ M " ds - j ] c(s)e Ar 'Hi 2 (s) ds (4j a ) 

— j I c(s)e~ " ds 
Jo 

+ ,/ f c(s)e Ars ds f f c(t)c(u)e- 1Vu+u) dt du 

J n J n J n 



H n (z) = j (' c(s)e~ M '' ds - j [ c(s)e ira H n _ l 2 (s)ds 
Jo Ju 

A (») = 

A,(«) = -j [ c(s)c Ars // 1 (s) ds 
Jo 

= f c(s)e* T ' ds f c(0« _Ar ' dt 

J{) Jo 

A 2 (z) = -j [ c(s)e sr °H 2 (s) ds 
J o 

= f c(s)e AT 'd8 f c(t)e~ ATt dt 

Jo Jo 

+ f c(s)c Ar " ds f c(t)e Ari dt [ [ c(u)c(v)e~ mu+v) du dv 

Jo Jo Jo Jo 

A.(f) = -j [' c(s)c AV °H n (s) ds. 



(41b) 
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Note that Ai is identical to the approximation of (26b). Writing as 
before 



A,,0) - A^O) = X„0), 
H n (z) - H n ^(z) = h n (z), 



we have 



(42) 



(43) 



44b) 



A H (Z) = EU2), 
fc=l 

H n (z) = £/>*(*). 
k=i 

The quantities \ n (z) and h n (z) are given as follows: 

X„( 2 ) = -j f c(s)e AVs h n (s) ds, n ^ 1. (44a) 

h n (z) = -j f c(s)e AT VI n ^ 2 (s) - 7/ H _ 2 2 (.s)] ds 

= -j f c(s)e Ar % ll . 1 (s)[H H - 1 (s) +H„-,( S )]ds,n ^ 2. 
•'o 

h(z) = H t (z) = j [~ c(s)e- AV ° ds. (44c) 

Jo 

Then under certain conditions described below, 

A(«) - E *»(«), (4.5a) 

n-l 
n-l 

We next obtain bounds on | \„(z) | and [ fe„(2) | . As stated in the 
last section, it is first necessary to impose additional restrictions on the 
problem. We assume that the coupling coefficient c{z) and the differ- 
ential loss Aa are such that a number K exists satisfying the following 
relation : 

f | c(s) | e Aa(z ~ s) ds ^ K for every z ^ 0. (46) 

We recall from (3) that Aa < 0. It will subsequently appear that 
convergence of the approximate solution can be guaranteed in general 
only for K ^ 0.455; further, the smaller K the more rapid the bounds 
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on | \„(z) | and | h n (z) | decrease as n increases. The restriction of (46) 
was again obtained partly by physical reasoning and partly by trial and 
error. Roughly speaking, for a given K it guarantees that the coupling 
coefficient c{z) is more or less uniformly distributed along z; we thus 
rule out cases where the coupling coefficient is zero over most of the 
guide and large over a very short section (e.g., the example of (31)-(34), 
where c(z) is a single 5-f unction ). Physically, such a condition says that 
for small K the spurious mode is dissipated much faster than it is coupled 
from the signal mode; the larger c(z) the larger must be | Aa \ in order 
to satisfy (46) for a given value of K. This will normally be the only 
case of practical interest in long random guides. 

Bounds on the first few h n (z) are readily obtained. From (44c) and 
(46), 



lh(z) | ^ f | c(s) I e~ Aa8 ds = c~ Aa: f 



c(s) 



Aa(z-s) 



ds 



< V 



- Aa: -K. 



Next, from (44b), (43) and (41a), 

h n (z) = -j2 f c( s )^ r X-i(s) f Zhis) +^ 



-i(s)' 



(47) 



ds. (48) 



Thus, 



h„(z) | ^ 2 f | c(s) | c Aa * | /»„_,(«) | [£" I h*( 
Jn \_k=\ 



s) 



+ 

Equation (49) yields for the first few h„(z) : 

I h(z) | ^ 2 f | c(s) | e la *e-' 21a8 ^ ds = e'^'K 2 
Jt\ 2 



h„-i(s) |~| , 
^~ J ds - 



(49) 



(50) 



SaU-„) , ^ ..-^ a - I' 3 



_/ \ I jo j-i j _» -jo: r- 

c(s) \ e ds ^ e -K 



ds 



h(z) | ^ 2 j" | c(s) | e Aas r^/OV^ 8 [~/y 4- 5. 

' '2/v 4 Tl 4- ?H ^' | c( S ) | C Aa, - s) rf s (51) 



= e 



< c "-2ZC 



1 + 



(53) 
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By an exactly similar process : 

| Uz) I ^ e-^-tfK 1 [l + |] [l + K 2 + K* + £] (52) 

I h(z) I ^ e- A -.2 3 K 9 [l + |] [l + K 2 + tf 4 + |] 

•|"l + K 2 + 2K 4 + 3/v 6 + 3K 8 + 3K 10 + 2K 12 + ^-1 . 

It is difficult to continue the above process and write out explicitly 
the n th term. However, by accepting a slightly poorer bound the analysis 
may be greatly simplified. We show that 

I K(z) | 52 e~ Aaz M"- 2 K n+l ; n ^ 2. (54) 

M is a constant to be determined, as a function of K. Assume that (54) 
is true for some value of n. Then, from (49) 



<t-i 

fc-2 trk+l 



I h n+1 (z) I ^ 2 f | c(«) | M"- 2 /r + V Aa '' X + E M k ~ 2 K 

JO _ k=2 

+ — — ] ds 

1 + K 2 Z M k K k + ^-A. (55) 

fc=0 ^ J 

• f | c(s) | e Aa(z - s) ds 

Jo 

^ e -^.2M"- 2 /C" +3 [l + K 2 g M fc K fc + M!_^'J . 
If (54) is to remain true for n —> n 4- 1, we have from (55) 

L A=0 ^ J 

But since the left-hand side of (56) is increased by allowing n —* °° and 
dropping the final term inside the brackets, the inequality of (56) will 
be satisfied if 

K \l + ^ 1<M. (57) 

L 1 - MKJ ~ 2 



(56) 
"" 2 ' 



APPROXIMATE SOLUTIONS 



1025 



M 



g.o 


















V 
V 

I 












1.8 


















I 
I 






















K + 


K 3 


M 

2 














1-MK 


1.6 
















> 




\ 

\ \ 


























N 












1.4 




















1 

1 












K I K ' < M 
'' ' 1-MK 2 

TO LEFT OF 

AND CURVES 
















1 


f 1 - 


Vl"(2K) 


2 


1.2 








II 


1 


K 

SPECIAL CASE 
C(Z) = C 
A/i = _ 
LARGE Z 








§ 


1.0 
















































0.8 


























































0.6 




M = 2.225K 

(POORER BOUND 

FOR K<0.3) 






















y 


















0.4 








\ 


// 
























/, 


// 






















0.2 




/ 


// 
s 


























y 





























/ 





























O 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

K 



Fig. 1 — M vs K 



A plot of M vs K taking (f)7) as an equalit}' is shown in Fig. 1; the 
inequality of ( .17 ) is satisfied to the left of this curve. If for a given A' 
we have chosen .1/ to satisfy the inequality of (57), then since (54) holds 
true for n = 2, (50), it is valid for all n by induction. For a given K we 
should choose the smallest value of M satisfying (57) in order to obtain 
the best bound. This smallest value of M is given by the solid curve of 
Fig. 1 (i.e., M < 1.554); the other branch, indicated by the dashed 
curve (i.e., J\I > 1.554), thus has no significance for our problem. We 
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note that convergence of the series solution of (45), and hence of the 
successive approximations of (41), is guaranteed only for ^ K ^ 0.455; 
for greater values of K the present analysis cannot guarantee con- 
vergence. 
Summarizing the above results: 

I h(z) | g e~ Aaz -K 

| h n (z) | ^ e-^-M"-' 1 K" +i ; n ^ 2. (58) 

^ K % 0.455. 

M is given as a function of K by the solid curve of Fig. 1. If A' is further 
restricted and if we are willing to degrade the bounds slightly, their 
form becomes simpler still. For example, if ^ K ^ 0.3 we may replace 
the bound of the solid curve on Fig. 1 by the slightly poorer dotted 
chord drawn from the origin. For this chord M = 2.225 K and the 
results of (58) become: 

I h(z) | ^ e" Saz -K 

I h„(z) | ^ e'^-K 3 (2.225 X 2 )"" 2 ; n ^ 2. (59) 

OgXi 0.3. 

Finally, by (45), (57) and (58) 

I H(z) I =g e~ la 



■* + £jr^]-.-[* +r ^]. 



(00) 



-*. M 



I H(z) | ^ e 

where M is again given as a function of K by the solid curve of Fig. 1. 
Having found bounds on h u (z), we may now find bounds on X„(z), 
our original objective. From (44a), 

I K(z) I ^ f I c(s) | e Aaw | /i„(s) | ds; n ^ 1. (61) 

Jo 

From (58): 

| Xx( 2 ) \ £K [ \ c(s) | da. 

I X„(z) | ^ M"- 2 K n+l f | c(s) | d«; n ^ 2. (G2) 

| A(z) | ^ M P | c ( 8 ) | ds. 
Ji Jo 

< K < 0.455. 
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M is again given by Fig. 1. Again, if K is further restricted, simpler but. 
slightly poorer results are obtained. For example: 

| \ x {z) | ^ K [ | c(s) | ds. 

I \ n (z) I ^ K\222r,K-)"~ 2 [ I c(s) I ds; n ^ 2. ( 63 > 

^ K ^ 0.3. 

Finally, the slightly better bounds of (51-53) may be used for the 
smaller values of n. 

We may again ask whether these bounds are the best that can be 
obtained. The answer is that we might be able to do a little better, but 
not much. Thus, consider the following special case: 

A/3 = 0, 

(64) 
c(z) = Co = pure real. 

From (46) we have 

K = -4- . (65) 

— Act 

The solution to the coupled line equations, (35), subject to the initial 
conditions of (36), for this case may be written in the following form. 1 

g,m = -' -Vi-wr- cxp Aa (1 + Vl _ (2Ky] 2 

2 Vl - (2 ^ 2 (00) 

For K < 0.5, all of the radicals in (66) are pure real. Under these con- 
ditions the first term of (66) has a smaller coefficient and a more rapidly 
decaying exponential factor than the second term. Therefore, for a large 
enough value of z the second term dominates, and we may write 



ft(.) » ' ± /' ~ f 2 */ cxp £ [1 - Vl - (2K)«fc (67a) 

2 \/ ] — (9k V 
A = -ln(?o(2) » -In 



2\/l - (2/C) 2 

1 + Vl - (2/Q g 

2 vr^w (67b) 

- ^ [1 - Vl - (2/v)%; 
K<0.5, -Aaz» vr ±^=. (67c) 
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The first term of (67b) is constant and the second increases linearly 
with z, so that we may write for large z 



large z. 



The bound of (62) becomes simply 



|A|^c„3. (69) 

Therefore, a comparison between the bound of the present analysis and 
the exact results for the special case of (64) for large z is obtained by 
plotting 

1 - \/j - (2/Q 2 
K 

on Fig. 1 and comparing this quantity with M . 

We see that for K < 0.36, the exact solution is indistinguishable from 
the bound of (69) on the plot of Fig. 1; consequently, little improve- 
ment may be obtained in these bounds without further restricting the 
problem. We also note that for K > 0.5, the above approximations made 
in the exact solution of (66) no longer apply. For K > 0.5 the character 
of the solution changes from monotonic to oscillatory; G Q (z) now has 
periodic zeros, at which In G (z) must approach infinity. Consequently, 
the series expansion for A in this case will diverge for K > 0.5. The 
present analysis guarantees convergence only for K < 0.455; while 
this might be a little smaller than necessary, the series solution may 
diverge for values of K not much larger. 

IV. DISCUSSION 

If K of (46) is very small compared to 1, jf« 1, the bounds of (62- 
63) on | X„(z) | converge very rapidly. Under these conditions it is 
tempting to assume that A (z) is satisfactorily approximated by the 
first term of the summation of (45); i.e., from (41b), 

A « Ax(aO = f c(s)e* r °ds f c(«)e" Art dt, (70) 

or one of the alternative forms given in (27). [Alternately we might 
wish to make a similar statement for G*(z) t as in Equation (24), when 

/ | c(s) \ds « 1.] This assumption has been made in all calculations of 

Jo 

transmission statistics that have so far been made. 
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Unfortunately there is,' at present, no satisfactory justification for 
this assumption. If | Xi(z) | turns out to be equal to its bound, as given 
in (62), and if K « 1, then, of course, we are guaranteed that the higher 
terms will have much smaller magnitudes than the principal term Xi(z). 
However, this situation is quite improbable, and occurs only in very 
specially selected cases. Thus if | \i(z) | is much smaller than its bound, 
as will be the usual case, we have no assurance that the magnitude of 
the next term | X2(z) | or higher terms may not be much greater than 
| \i(z) | . However, no instance is known in which | X 2 (z) | is not small 
compared to \\i(z) \ , for K « 1. 

We do not know whether or not the perturbation solution of (70) 
provides a useful approximation for all cases of interest (i.e., for all 
cases where c(z) satisfies (46) for some small value of K, e.g., K = 0.1 ). 
Even if this approximation fails in some cases, we may still hope that it 
holds true in most cases, so that (70) will yield the correct statistical 
properties of the loss and phase when the coupling coefficient is a sta- 
tionary random process, perhaps Gaussian, with a sufficiently small rms 
value, at least for the simpler statistics of interest. Although this is be- 
lieved to be true by a number of people, there is nothing in the present 
paper that bears on this question (and no other information known to 
the author). It would be most desirable to obtain further information 
on the way in which Ai(z) of (70) approximates the true solution A(z); 
e.g., does Ai(z) approximate the fine structure of A(z) as well as its 
average value as A/3 (which varies with the frequency of the applied 
wave) varies. 
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